Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis     

Ramin Pakzad  Shanyong Wang  Scott W. Sloan 《国际地质力学数值与分析法杂志》2020,44(10):1472-1494

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.  相似文献   

Solution to the rotation of the elastic earth by method of rigid dynamics     

Yoshio Kubo 《Celestial Mechanics and Dynamical Astronomy》1990,50(2):165-187

Hamiltonian mechanics is applied to the problem of the rotation of the elastic Earth. We first show the process for the formulation of the Hamiltonian for rotation of a deformable body and the derivation of the equations of motion from it. Then, based on a simple model of deformation, the solution is given for the period of Euler motion, UT1 and the nutation of the elastic Earth. In particular it is shown that the elasticity of the Earth acts on the nutation so as to decrease the Oppolzer terms of the nutation of the rigid Earth by about 30 per cent. The solution is in good agreement with results which have been obtained by other, different approaches.  相似文献   

有限元法与伪谱法混合求解弹性波动方程   总被引：6，自引：0，他引：6  

马德堂  朱光明 《地球科学与环境学报》2004,26(1):61-64

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.  相似文献   

基于Poulos弹性理论的被动桩改进算法     

陈福全  黄伟达 《岩土力学》2008,29(4):905-910

传统的Poulos弹性理论仅适合于均质土中土体侧向位移时桩的性状分析,无法考虑土的层状特性。通过引入层状地基中作用一水平集中力的广义Mindlin解和地面作用有竖向荷载时的应力和位移通解,对Poulos方法进行了改进,使之扩展到多层土中,还用于研究堆载条件下的被动桩变形和受力响应。算例分析表明,改进弹性理论要比Poulos方法更为严密、合理,提高了计算精度,应用范围也更广。  相似文献   

Reduction of hydroelastic responses of a very-long floating structure by a floating oscillating-water-column breakwater system     

D.C. Hong  S.Y. Hong  S.W. Hong 《Ocean Engineering》2006,33(5-6):610-634

The hydroelastic responses of a very-long floating structure (VLFS) placed behind a reverse T-shape freely floating breakwater with a built-in oscillating water column (OWC) chamber are analyzed in two dimensions. The Bernoulli–Euler beam equation is coupled with the equations of rigid and elastic motions of the breakwater and the VLFS. The interaction of waves between the floating rigid breakwater and the elastic VLFS is formulated in a consistent manner. It has been shown numerically that the structural deflections of the VLFS can be reduced significantly by a suitably designed reverse T-shape floating breakwater.  相似文献   

作用于平台桩腿的海冰屈曲破坏冰压力随机过程模型及其参数确定   总被引：2，自引：1，他引：2  

段忠东  欧进萍 《海洋工程》1996,(4)

冰对结构的作用过程是典型的随机过程。本文在冰压力过程为平稳过程的假设下，从渤海海冰对平台桩腿作用的大量实测冰压力数据中，选取了２１条冰屈曲破坏时的冰压力时程曲线，对它们进行谱分析后，得到了单点冰屈曲破坏的压力随机过程的谱密度，并确定了谱参数及其跟环境要素的关系，依据文献［１］中冰压力沿圆柱面的空间分布，建立了绕桩腿的冰压力随机场模型，并得到了作用于桩腿的总冰力随机过程及其谱密度。本文的研究成果为平台结构冰激随机振动和疲劳累积损伤分析提供了荷载基础  相似文献   

Impact Buckling and Post-Buckling of Elasto-Viscoplastic Cylindrical Shell Under Torsion     

Wang Deyu Chen Tieyun Associate Professor 《中国海洋工程》1997,(1)

The dynamic torsional buckling and post-buckling of elasto-plastic circular cylindrical shell isstudied.By the usage of the Bodner-Partom constitutive relation,the present visco-plastic problem is treat-ed as an orthotropic elastic problem.Based on the nonlinear large deflection shell therory,the governingdynamic buckling equation is solved by the Runge-Kutta method,and the critical buckling stress is deter-mined by the B-R criterion.The present paper is mainly focused on the imperfection sensitivity,the strainrate sensitivity,and the dynamic post-buckling characters for the shell loaded dynamically with a constantshear strain rate.Numerical results show that the critical shear stress and imperfection sensitivity will in-crease with the increase of strain rate.It is found that the elastic buckling stress is more sensitive toimperfection than the visco-plastic buckling.The present research also reveals some dynamic post-bucklingcharacters for the cylindrical shell,and it is shown that the shell will vibra  相似文献   

Analytic Solution of Interaction Between Surface Waves and Floating Elastic Plates with Elastic End-Restraints     

赵存宝张嘉钟  邢海燕黄文虎 《中国海洋工程》2007,21(1):23-38

Based on the dynamic theories of water waves and Mindlin plates,the analytic solution of interaction between surface waves and two-dimensional floating elastic plates with edge-restraint is constructed by use of the Wiener-Hopf technique.Firstly,without regard for elastic edge restraint,the wave-induced responses of elastic floating plate analyzed by the present method are in good agreement with the results from literature and experimental results.Therefore,it can be shown that the present method is valid.Secondly,three end-restraint cases(i.e.,the left-end elastic restraints,the both-end elastic restraints,and the right-end elastic restraints)are proposed to reduce the vibration of floating plates,in which the spring is used to connect the sea bottom and the floating plate's left(or right)edge.The relations between the spring stiffness and the parameters of wave-induced responses of floating plates are discussed.Moreover,the effective method to reduce the vibration of floating elastic plates can be obtained through comparison.  相似文献   

波浪作用下潜标拖缆动力学分析     

朱刚  杜月中 《海洋工程》2007,25(4):15-20

对三维潜标拖缆在规则波作用下的动力响应进行了研究,首先基于集中质量法给出了缆绳在规则波作用下的控制方程,然后提出了两种不同类型的自由端边界条件:质点型和艇型自由端边界条件,前者把拖体视为作三个自由度线运动的质点,而后者把拖体视为作六自由度运动的艇体,最后用有限差分法进行了数值研究。数值计算结果表明潜标拖缆在规则波的作用下,不仅作近似的简谐振动,还会在横向、垂向产生一个偏移量,该偏移量的大小随浪向角的变化而变化,另外受遭遇频率的影响,振动的幅值也随浪向角的变化而呈现不同的特征。  相似文献   

海冰板的环形屈曲的压缩破坏   总被引：2，自引：0，他引：2  

唐立强  陈会军等 《海洋工程》2001,19(4):29-33

在轴对称条件下，给出海冰板的环向屈曲面临载荷，从细观力学的观点分析了微裂纹开裂的条件，扩展方向并给出绕海洋平台柱腿冰的损伤区的尺度。所给出的临界半径是一个与冰材料力学特征，冰的微结构参数一海洋平台柱腿半径相关的特征参数。对现场所观察和测量的结构给出了解释和定量分析。  相似文献